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Let P=\p n \ be a real sequence and define the transformed sequence Q (P) as follows: 
Consider the segment, p m , p m +\, p m +2, Pm+z of P and determine a m , b m , and c m so that 
m.Q J x\a m bZ J r^n — Pn\, n = m, ra+l, m-f-2, m + 3 is minimized. The rath term of f (P) 



n 
is taken to be c, 

k k I 

Po-r-S^-^i P° + (- 1) n Xl °wVX! M 



The effect of the transformation C ° n sequences of the type 

a\ n . 



and po-\ is considered. In each case £ * s 

.-l i=i i=\ n 

shown to be very effective in accelerating convergence or decelerating divergence. For 
example, if the second sequence behaves as an k as n— >co ; then the transformed sequence 
behaves as an k ~ 3 . A similar transformation T> is defined by approximal ing in the Tcheby- 
cheff sense a segment of P by ab n cos (0-\-n<p)-\-c. The effect of T) is studied for sequences 
of the above type and also for po + 2t=i«i^? cos (0. + w0t)- These sequence transformations 
are similar in nature to Aitken's 5 2 -process and its generalization. A comparison of the 
two types of transformations is made. Several examples are given to illustrate the effect 
of C and T) on various sequences. 



1. Introduction 

Let P={p n \n = l,2, . . .} be a real sequence and 
let T be a sequence transformation with 



Q=T(P) = {t n (p hP2 , 



p n )\n=l, 2, 



A large variety of such transformations have been 
studied, many in the modified form of summability 
methods for infinite series, which are designed to ac- 
celerate or induce the convergence of the sequence P. 
Such transformations may be based on approximat- 
ing a segment of P by a function of a given form. As- 
sume that p n behaves approximately asJ(A f n) where 
A represents some parameters to be determined. 
Consider a segment {p i \i=-l,2, . . . , m] of P. The 
parameters A m may be determined so as to minimize 
the error of the approximation in some sense. Since 
f(A m ,n) is a known function of n we may set 

tm(Pi,P2, . . ., p m )=IAmf(A m ,ri). 

The transformations of this paper result from as- 
suming p n to behave as ab n +c or ab n cos (d+n<t>) + c. 
The transformations are denoted by Q and *T), re- 
spectively. The approximation to give the mth 
term of C(P) an d T>(P) is based on four or six, respec- 
tively, terms from P. The parameters are determined 
so as to minimize 



Tnax\a m b£+c m —pn 



n=m, m+1, m+2, m+3 



or 



max\a n bm cos(d m +n(p m ) +c m —p n \, n=m, . . ., m + 5 



1 Present address, General Motors Corp., Warren, Mich. This work was 
done under a National Research Council-National Bureau of Standards Post- 
doctoral Research Associateship. 



If 6 OT <0, then Lim j(A m ,n) is just c m ; if b m ^>l then 

W->oo 

this limit does not exist, but the value c m is still taken, 
c m being the antilimit. 

Explicit formulas for c m in terms of the p t have been 
obtained in [5]. 2 These formulas are not linear in the 
Pi and the transformations exhibit some marked dif- 
ferences from the classical linear transformations, in 
many cases Q and T> arc much more powerful and in 
other cases P may converge quite rapidly and the 
transformed sequence may be identically infinite. Al- 
though C and *T> are not regular they do have a prop- 
erty of joint convergence, namely that if both P and 
the transformed sequences converge then they have 
the same limit. 

It is not surprising that Q and T) are effective for 
sequences of exponential type 



'-{ 



Po+Z)«iX?|»=l 



-■•■}• 



(1) 



If |Xj|>|X i+1 | then C(P) converges or diverges as 
(X 2 ) w instead of (XO*. T)(P) behaves as (X 3 ) w or 
(a!Ai)*- T) is effective for sequences of the more 
general type 



: /^+Z>*X? 



COS (0* + 710 <)| ft 



= 1,2, . . .}• 



(2) 



and C is useful for such sequences in some cases. 

These transformations are also very powerful for 
alternating rational sequences: 



^ i=l i=l J 



(3) 



1 Figures in brackets indicate the literature references at the end of this paper. 
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It is shown that £ accelerates the convergence by a 
factor of 1/n 3 . That T) accelerates the convergence 
even more can be seen from the examples although 
the actual factor of acceleration has not been deter- 
mined. An example shows that £ an( i ^ are very 
effective for a sequence of the form 



.{i=23! ,„-!,*,...}. 



(4) 



The basic principle of these transformations is the 
same as that of Aitken's <5 2 -process and its generali- 
zations. £ i s the direct analog of the <5 2 -process. 
Shanks [8] has given the most complete account of 
this theory. £ an d ^ have many properties in com- 
mon with these transformations. The analogs of £ 
and *D in Shanks' work have the same power as £ 
and T)for sequences of type (1) and (2). For se- 
quences of type (3), C an d ^D are definitely more pow- 
erful. Aitken's 5 2 -process accelerates the conver- 
gence of (3) by a factor of only 1/n 2 . £ an d *2) are 
much more effective for the sequence (4). There are 
sequences for which the <5 2 -process and its generali- 
zations are more effective, but as a general guide £ 
and T) are as effective or more effective than their 
analogs in Shanks' paper. 

A comparison of this paper with Shanks' will indi- 
cate that £ and T) have many properties which have 
not been developed. The main purpose here has 
been to establish these transformations as tools of 
numerical analysis and not of analysis. When read- 
ing Shanks' paper one has the feeling that much 
remains to be done in the study of these transforma- 
tions and that a complete understanding of their 
behavior for sequences of real and complex functions 
will require a penetrating analysis. 

It is typical but somewhat disappointing that these 
transformations may be most effective for easy 
problems and least effective for difficult problems. 
If one is solving three simultaneous linear equations 
by Gauss-Seidel iteration and the characteristic 
values determining the rate of convergence are .95, 
.5, and .1, then these transformations increase the 
rate of convergence dramatically. But if one has 
100 equations with 10 characteristic values between 
.99 and .995, then the elimination of even the five 
largest ones is not very significant. On the other 
hand if the largest characteristic value is .995 and the 
others are, say, less than .9 the transformations are 
very effective in accelerating the convergence. 

2. Derivation of the Transformations 

2.1. Preliminaries 

Sequences are denoted by P, Q, . . . and their ele- 
ments by p n , q n , . . . . A sequence transformation 
is denoted by a script letter as £ and T). Braces, 
{ }, denote a sequence or set and {x\ . . .} is read 
"the set of x such that . . ." We denote by 0(x) 



and o (x) two functions such that 



and, for x sufficiently large, 



0(x) 



< constant. 



Let P be a sequence of real numbers. A trans- 
formed sequence Q={g n \n=l, 2, . . .} of P, is ob- 
tained as follows: Take a segment, {pi\i=m, . . ., 
m+&} of P, and approximate the values p n as a 
function of n by ab n +c or ab n cos (6 +n4>)+c. The 
corresponding c value is then assigned to q m . 

For these particular transformations to be effective, 
the sequence P must, in some sense, behave ex- 
ponentially. As an example of such sequences 
consider the real ^-vectors defined by 



Vn+l- 



~-Av n) 



where A is a real kXk matrix. Let v* be the solution 
of 

Then v n+l — v*=A n (v — o*). The sequence of vectors 
defined will converge to v* it A n (v Q —v*) tends to zero 
as n tends to infinity. If A is a normal matrix with 
characteristic vectors u t and distinct real character- 
istic values X f , then 

If the A f r are complex then 

A n (v -v*)=J2^i\\i\ n cos {e i +n4> l )n i . 

2.2. Tchebycheff Approximations 

We wish to approximate a segment of P by ab n -\-c. 
The use of least squares approximation appears to 
be impossible due to the difficulty of the nonlinear 
equations involved. In [5] the theory of approxima- 
tion by ab x +c in the Tchebycheff sense is developed 
in detail. The following result is given. Let p m , 
p m +i, Pm+2, and p m+3 be four consecutive values of P 
to be approximated. Then c m of the best approxima- 
tion, a m (b m ) n +c m , is given by 

c ^ (Pm+Pm+l)(Pm+2+Pm+s) — (Pm+l + Pm+2) 2 , g v 
2(Pm— Pm+1— Pm+2 + Pm+z) 

Let p m , . . . , p m +5 be consecutive values from P 
to be approximated by a& w cos(0+7i0)-|-d. The theory 
of Tchebycheff approximation by ab x cos(d J rx<j))+d is 
not complete and an explicit formula for d m of the 
best approximation has not been rigorously estab- 
lished. However it is conjectured that the following 
formula is valid. Set s i = i / 2 (p m+ i J rp m+ i +1 ), z'=0, 1, 
. . ., 4. Then 



2S1S2S3-}- §4(50^2 ^1/ S()$3 S 2 



2s 2 (s +2s 1 —Ss 2 J r2s 3 +s 4 )—2(s 1 —s 3 ) 2 —4:(s l s i +s s 3 ) + 2s o s 4i 



(6) 
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This formula is based on the assumption that the 
best Tchebycheff approximation is characterized by 
the alternation five times of the error function. 
See [6] for a discussion of the characterization of 
best Tchebycheff approximations. 

Note that (6) may be written in a much simpler 
form. Let s , i =(s i —s 2 ) y i=0,l,3,4; then 



d m — s 2 



s'tisiy+s'o&y 



5 o^4 \ s i s 3 ) 2 (Si<s 4 ~r s^Sq) 



2.3. The Transformations Q and T) 

Let P be given, then a new sequence Q is found as 
follows : For each m such that the denominator of (5) 
is nonzero q m =c m ; if the denominator of (5) is zero 
and the numerator does not vanish then we assign 
q_m= °° ; if both the numerator and denominator are 
zero q m =Pm> This transformation of P into Q may 
be written in operator form as 

Q=C{P). 

A transformation T) is defined in a similar manner 
by (6). 

The repeated application k times of C and T> is 
denoted by C(P) and T> k (P). 

3. Analysis of the Transformations 
3.1. Algebraic Properties of Q 

The study of the properties of these transforma- 
tions begins with a simple algebraic property of C- 

If P={p n } then aP-\-b is defined as {ap n J \-b\n= 1,2, 
. . .}. The following result maybe established by 
direct computation from (5). 

Theorem 1: Q (aP+b)=aC (P) + b. 
Itis not true that C(Pi+P 2 )=C(Pi)+C(P2). 
3.2. Transformation of Exponential Sequences by Q 

Let P be a real sequence of the form 

\ Pn=Po+lt l <x i \"\n=l,2 f . . . V 

with |Xf|>|X z+ i|. If I X 2 J ^> 1 , Pdoes not converge but 
diverges from p . In such a case p is the antilimit 
of P. This sequence is of an exponential type and 
C should be effective in increasing the rate of con- 
vergence or decreasing the rate of divergence. 
Let P'=P-p ; then 

C(P)=C(P')+Po. 



The general term of C(P f ) niay be explicitly com- 
puted to be 

i<i 

a 1 (l-.X 1 -X?+X?)+S«i(X < Ai) n (l-^-M+X?) 



The denominator of this expression may have only a 
finite number of zeros as a function of n. We may 
write the denominator as 

ai(l-X 1 -X?+X?)+o(l) 

since (X^/XiKl for j>l. 

The largest term in the numerator is X*. Others 
among the larger terms are 

x?, (WXi)". 

Thus C is seen to eliminate the largest exponential 
term from P'. 

The repeated application of C will eliminate the 
largest remaining term. The largest term in C 2 (P') 
is AJ. 

The above analysis has established: 

k 

Theorem 2: Let P={;po+2>iM}, |Xi|>|X i+1 |; 

i=i 
then 



C(P) = iPo 



1 + 0(1) 



Mo=l 



Thus C is seen to have a desirable effect on se- 
quences of this type, which was to be expected. In 
the next section it is seen that C is effective for 
some sequences of a completely different nature. 

3.3. Transformations of Alternating Rational 
Sequences by Q 

Let P be a sequence of the form 

^ t = / i=0 J 

We shall consider C(P)* 
First take the special case 

P={po+(-l) n an k }. 

The nth term of C(P) is 

(n k -(n + l) k )((n + 2) k -(n4-3) k )-((n-{-l) k -(n-\-2) k )z 



Voi-a- 



2(-l) n (ri k +(n+l) k -(n + 2) k -(n + ?>) k ) 
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After some manipulation this term is seen to be of the 
form 

an 2lc -*k 2 (k— 1) . p( k-±\ 



or 



8(— l) n kn k 

,0 + .(-DV[^+O([ 4 )]. 



Now consider the general situation. p n may be 
written as 



Pn=Po+o&n k + y Z2a[n? 



k=ki—k 2 



for n sufficiently large. It is clear that the denomi- 
nator of c m is 8{-l) n kn k - l + 0{n k - 2 ) since it is linear 
in the terms of P. The denominator of c m has only 
a finite number of zeros as a function of m. The 
terms in the numerator are of the form 

«W[n w - (n+ 1)*-'][(ti+2)*->- (n+3)*-1 

+ ^a;[^^-(^+l)^][(7i+2) A - i -(7i+3) fc - i ] 

By direct calculation it is seen that this term is 

0(n 2k - *-*-*). 

Thus we have established: 

Theorem 3 : Let P= {p«+ (- 1) "S*-aa<ft7Si-oM f } ; 
^en /or m sufficiently large the mth term of C 
(P) is 

OO 

2=1 

where k=k 1 —k 2 , a =a kl k(k—l)/Sb k2 . 

Examples of sequences which are approximately of 
this type are given by the partial sums of alternating 
infinite series of rational functions of n. One can 
consider such sequences directly and establish an 
analog of theorem 3. 

3.4. Regularity and Joint Convergence 

A transformation T is said to be regular if the 
existence of Lim p n implies Lim g w =Lim p n where 

n-$ oo n-> oo w— > oo 

Q=T(P). Since £ is a transformation designed to 
accelerate or induce convergence one would like for 
C to be regular. That C is not regular is seen by 
the following example: 

P={Pn\p2 k =l + lItc, P2ic+i=l + l/k}- 

Every other term of C(P) is infinite, yet P con- 
verges. Although C is not a regular transformation 
it does have the property of joint convergence. A 
transformation T is said to have the property of 
joint convergence if the convergence of P and T(P) 
imply that P and T{P) have the same limit. The 
proof of the following theorem will be given in section 
.5.1 (see theorem 10). 



Theorem 4 : If both P and C(P) = Q converge then 
Lim p n = Lim q m - 

W-> oo ?/i-> oo 

3.5. Properties of the Transformation T) 

The following theorem gives algebraic properties 
of T) similar to those of C- The proof is by direct 
evaluation. 

Theorem 5: D(aP+b)=a r L>(P)+b. 

Consider the real sequence 

P=/]Ca*A? cos {B i +n<t> l ) \ 

with |Xi|>|X f+1 |. Expression (6) may be directly 
evaluated as the ratio of the following two terms : 

S a i a j a l (\ i \ j \ l /\ 2 1 ) n 
!*-;H-li-/|>o 



and 



Xcos [6 i +e j +d l +n(ct> i +cl> j +<f> l )]H ijl 



2 = 1, .7=1 



where H tj i and G tj are functions of the X's, 0's, and 
0's that are independent of n. H in and G tJ are 
trigonometric functions with coefficients depending 
on powers of the X's. 

An examination of the above expression and a 
simplification of notation leads to the following 
theorem. 

Theorem 6. Let P={po+^ k =1 ai\ n i cos (O t +n<l>i)}, 
|X<|>|X<+i|. Then the nth term ofT){P) is 



Po~ 



X2 S &<M? COS (ai+nfii) _ -. 

z=0 Mo— J- 



1 + 0(1) 



H\<h t>l 



except for possibly a finite number of terms. 
One may also establish the following result : 



Theorem 7: Let P={pv+l>Za i \ n i }, |X*|>|X <+1 1. 
Then the nth term ofD(P) is 



Po+wmb^ 

2 = 

1 + 0(1) 



Mo=l 

l/*,KM>i 

X' = max (X 3 , Xl/Xi) 



except for possibly a finite number of terms. 

No attempt has been made to analyze the effect 
of T) on alternating rational sequences; however the 
examples indicate that T> is quite effective for in- 
creasing the rate of convergence. 

It is a simple matter to construct sequences for 
which T) is not regular. It is probable that T) has 
the joint convergence property, but this has not been 
established. 



:3 



4. Comparison With Aitken's 5 2 -process and 
its Generalization 

The transformations C and 1) are similar to a 
family of sequenee transformations first studied by 
Aitken [1,2]. These transformations have been stud- 
ied since by Lubkin [4], Samuelson [7], and Shanks 
[8]. The most general and complete treatment is 
by Shanks. A brief outline of these transformations 
will be given. 

Let P be given and set 

Ap n =p nH — p n . 



From P new sequences Q k = 
are formed by 



{q k , n \n=k, k+1, 



9k, n 





Pn-k • 


. • Pn- 


I Pn 




APn-k • 


• • kPn- 


&Pn 




&Pn-k-l • 


. . Ay;, 


bPn+1 




bPn-l • 




&Pn+k+l 




1 


. . 1 


1 




&Pn k • 


. • &Pn- 


1 kp n 




kpn-k-1 • 


. • ^Pn 


kPn+1 




&V n -\ . 




&Pn+k-l 



This defines a set of sequence transformations £ k) by 

Q k =£ k (P). 

Shanks studies in some detail the transformations 

£k,£iA = {pi,q.i2,g.W, • • • , <&.£+i> ■ • • }> and 

£ d ={q k , k \k = 0,l,2 } . . . }. (iin } denotes the nth term of 
£ k (P). qin is given explicitly by 



Pln = 



Pn + lPn-l — PJ 



Pn+l+Pn 



-2p n 



which defines Aitken's <5 2 -process. 

These transformations are derived heuristically as 
follows: Given a segment of P, p n -k, • ■ ■ , Pn+k, 
assume that this segment is exactly of the form 

k 

2«jW+c Then define q kjn to be this c. The basic 

i = \ 

principle is the same as for C and 1), but the 
determination of the approximation is different. 
If S={s n \Sn=?(Pn+Pn+i)} then it is seen that 

C(P)=e l (S). 



Likewise, a straightforward calculation shows 
'D(P)=e 2 (S). 



The reason for this relationship is simple 
been found by solving eq (8) for c: 



8k h 



y £^a i bf+c=p n 



m = n - 



,n + k. (8) 



The fact that best Tchebycheff approximations are 
characterized by the alternation £+1 times of the 
error function leads to the following equations: 



^aM+c=p m +(— l) m e, ?n = n, . 



n+2k+2 



where e is the error of the best approximation. If 
one averages the mth and (m + l)st equations, then 



± ai{l + bi) bT 



"C=2 (Pm+Pm+l), 



m=n, 



n+2k+h 



These equations are of the same form as (8) except P 
has been replaced by S. 

Aitken considers the application of £ x to sequences 
of the form P={p t ) J r^2 k i = \a i X7} and establishes: 

k 

Theorem 8 (Aitken): Let I > ={po + ^ l ^iK} ) 
|\<|>|X <+ i|; then 

£i(P) 



Po- 



i=0 






1 i>l. 



Thus for sequences of t Ins type £\ and Q are equally 
effective in order of magnitude. The coefficients of 
(X 2 ) w for £i and Q are, respectively 



fti<3 2 (Xi— X 2 ) 2 
(1-XO 2 



«ig2(l + X 2 )(X 1 -\ 2 ) 2 

(i-xo 2 



£i appears to be slightly superior if X 2 >0. Examples 
indicate that either <5\ or C may be better but in any 
case the difference in the rates of convergence is 
small. 

For alternating rational sequences the two trans- 
formations are no longer similar in effect. Shanks 
has established: 

Theorem 9 (Shanks): Let 

P=(po+(-l)W/swV 

L *=o / *=« J 

dki^O, b k ^0 and k=ki—k 2 . Then the nth term oj 
ft(P) is 

a h k(-l) n n k - 2 



Po+ 



46* 



-0(n k -*) 



562759—60- 
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with possibly a finite number oj exceptions . 

Thus £ x accelerates the convergence by a factor 
of l/n 2 whereas C accelerates it by a factor of 1/n 3 . 

The comparative effectiveness of 8 2 and ^ mux 
be judged from the examples given later. 

Theorem 10 shows that <5\ has the property of 
joint convergence [4]. 

Theorem 10 (Lubkin): Ij both P and 8 l (P) = Q 
converge, then 

Lim p n =lAva q n 

n^> oo Tt — > oo 

If P converges then so does S= { }i(p n -hPn+i) } and 

Lim 2> w =Lim s n 
7fc— >oo n-^co 

Thus theorem 4 follows directly from theorem 10. 

5. Examples 

Some examples are given which illustrate various 
characteristics of C and ^- In. some cases com- 
parisons are made with <5\ and 8 2 . For more 
general comparisons see [4] where 8 X is compared 
with several classical linear sequence transforma- 
tions for a large number of examples. 

The first two sets of examples illustrate the general 
behavior of C and ^ for various exponential and 
alternating rational sequences. Some comparisons 



are made with £ x and 8 2 here. The next set of 
examples shows the effect of increasing the largest 
A in exponential sequences. The following two sets 
show the effects of lower order terms and the loss of 
significant figures. Detailed comments are given 
with each set of examples. 

All of these examples were computed using double 
precision arithmetic with about 20 significant decimal 
digits. 

5.1 . Effect of C an d *2) on Exponential Sequences 
The following five examples are considered: 

p 1 ={2 n +(i.2y+(.9y+cs) n }, 

P 2 ={(.98) n +(.95y+(.5) n +(.2) n }, 
Ps={(.98)*+(.95)»}, 
P 4 ={(.9y+(.Sr cos (.In)}, 
P 5 ={(.9) w cos (1.5n) + (.5) w cos (.25n)}. 

Table 1 gives results of C, C\, C\ C\ and H0> for 
33 terms of P x . I\ diverges quite rapidly and Q (I\) 
diverges at a much slower rate. The factors in 
C(Pi) larger than .5 are 1.2, .9, .8, .72, and .54. 
C 2 (I\) converges as can be seen from the list of 
factors in C 2 U\) larger than .5. They are .9, .8, 
.72, .68, .6, and .54. The four largest factors of 
Ca (P 2 ) are .8, .72, .71, and .68 are the largest factor 
of CKPi) is .72. 



TA.BLE 1. 



1\ 


C(Pi) 


C<Pi) 


CHPi) 


CKPi) 


*D(Pi) 


4.90000000000000 


2.9151478 


2. 7270889 


2. 7433842 


2. 7029023 


2. 6785648 


(i. 89000000000000 


2.7741418 


2. 7426402 


2. 8808482 


3. 0727475 


2. 7348679 


10. 9()90()()()()()( )()()() 


2. 7169114 


2. 6757849 


3.2921161 


69. 180737 


2. 7088105 


1!). 1393000000000 


2. 7430500 


2. 5368629 


4. 8080981 


72. 818042 


2. 6047836 


35. 4064900000000 


2. 8532800 


2. 3466523 


140. 43599 


-.54172894 


2. 4392636 


07. 7795(i90000000 


3. 0509855 


2. 1287500 


-2.4381772 


-.04037382 


2. 2345869 


132.271192900000 


3.3426610 


1. 9033466 


-.76363656 


.04995813 


2. 0125094 


260. 898056330000 


3.7381380 


1. 6847333 


-.28755651 


. 06710598 


1. 7902083 


517.681418569000 


4. 250855S 


1. 4814607 


-.09493555 


. 06620907 


1. 5791104 


1030. 64778904490 


4. 8982545 


1. 2976421 


-.00731646 


. 06136882 


1. 3855556 


2055. 82979364889 


5. 7023160 


1. 1344196 


. 03331326 


. 05701478 


1. 2121566 


4105. 26724946147 


6. 6902(194 


.99116645 


. 05084351 


. 05513287 


1. 0591391 


8203. 00848270213 


7. 8954805 


. 86633340 


. 05650783 


. 05808964 


. 92537882 


16397.1119330351 


9.3585517 


. 75798032 


. 05602878 


. 07372670 


. 80910184 


32783. 6480970788 


11.128667 


. 66412067 


. 05246752 


. 17953296 


. 70831457 


65554. 7018754061 


13. 265229 


. 58282968 


. 04750712 


-. 11327597 


. 62104494 


131094. 375400883 


15. 839833 


. 512382X9 


. 04208180 


-.03107621 


. 54546427 


262170. 791442315 


18. 938651 


. 45125306 


. 03670509 


-.01478933 


. 47993879 


524320. 097496628 


22. 665309 


.39811524 


. 03164886 


-.00839710 


. 42304196 


1048614. 47070579 


27. 144335 


. 35183192 


. 02704444 


-.00521395 


. 37354352 


2097198. 12376227 


32. 525314 


.31143263 


.02294119 


-.00341676 


. 33037177 


4194359.31199968 


38. 987872 


.27609212 


. 01934122 


-.00231783 


. 29269020 


8388674. 34190501 


46. 747669 


. 24510968 


.01622020 


-.00160176 


. 25998290 


16777295. 5813360 


56. 063584 


. 21789028 


. 01353993 


-.00114340 


. 23010570 


33554527. 4717843 


67. 246349 


. 19392818 


.01125615 


-.00081105 


. 20587142 


67108978. 5430931 


80. 668906 


. 17279283 


. 00932342 


-.00056122 


. 18756025 


134217865.431120 


96. 778842 


. 1541 16S7 


.00769611 


-.00044434 


. 14815635 


268435620. 898931 


116.11330 


. 13758611 


. 00633423 


-.00029985 


. 12460304 


536871109.862244 








. 00001497 




1073742061. 41994 












2147483932. 89072 












4294967637. 85702 












8589935002. 21781 
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It is seen that the first values of C k (Pi) behave 
erratically as k increases. C*(Pi) > s a very erratic 
function of n for small n. This is the typical behavior 

and indicates that care must be exercised in using 
repeated applications of C- 

The largest factor in T)(l\) is .9 and hence '/)(/',) 
converges. 

Only selected values from the other examples are 
given. In table 2 it is seen that C is n °t very effec- 
tive in increasing the rate of convergence of P 2 . 

Table 2 
P 2 ={(.98)»+(.95)»+(.5)«+(.2)»j 



11 


Pi 


CUV 


CKPi) 


r D(P i ) 


1 

10 
20 


2. 6300 
1. 4168 

1.(1261 


1.6876 
. 32965 
. L8278 


3. 0466 
. 20621 


1. 2528 
.21299 
. 02770 



This is due to the closeness of the two largest factors. 
In the next table it is seen that even with the lower 
order term deleted the convergence is quite slow 
even for C 4 (P 3 ). 

Table 3 

P 3 ={(.98)»+(.95)«] 



n 


Pz 


C(Pi) 


Pi/' > 


c /■ 


C P 


1><P«) 


1 


1. 93 


. 29914 


-.07157 


-.08733 


-.07864 


.0 


25 


. 88085 


. L5891 


.04624 


-4.9028 


. 36)843 


.0 


40 


. 57421 


. 10160 


-.02795 


. 02408 


.00322 


.0 


00 


.34362 


.05175 


-.01084 


.00256 


. 00004 


.0 



The next table shows that *D is fairly effective on 
P 4 . It is to be noted however that convergence is no 
longer monotonic even in the large range 4 considered. 
The denominator has a zero between the 67th and 
68th terms of r D(I\). This same phenomenon is 
present in T) 2 (P A ). 

Tabi e 1 
p 4 -{(.9)»+(.8)«cos(.ln)) 



n 


Pa 


C(Pi) 


1>(P0 


1P(Pd 


1 


1. 6960 


. 02354 


. 08971 


. 05260 


8 


. 54735 


.07107 


.11293 


-. 47672 


10 


. 40669 


. 06085 


-1.8972 


-.00681 


20 


. 11678 


.01132 


-.00282 


-.00287 


30 


.04117 


-. 00036 


-. 00051 


-.00048 


35 


. 02465 


-.00042 


-.00151 


. 00005 


50 


.00516 


-. 00002 


3X10-6 


2X10-7 


60 


. 00180 


4X10-8 


5X10-7 


3X10-7 


67 


. 00086 


3X10-8 


3X10-6 


6X10-7 



The final table shows that T) is very effective for 
P b , while C is n ot. 

Table 5 

Ps = { (.9) «cos (1 ,5n) + (.5) "cos (.25?/) ) 



n 


Ps 


OPs) 


1>(P 5 ) 


1)'(P5) 


1 
5 

10 
15 
20 
25 


. 54812 
. 21454 
-.26567 
-. 17983 
. 01875 
. 07037 


-.11751 
-.21059 
3. 2919 
. 07227 
. 06441 
-. 04624 


. 10538 
.00017 
-.00022 
-5X10-6 
1X10-7 
7X1O-0 


-.00135 
-1X10-8 

-6X10" H 
5X10-'i 
6X10-" 







5.2. Effect of C an d T) on Alternating Sequences 

Consider the sequences of partial sums of the follow- 
ing series: 

CO (_))> \\ 

- 4 §2i-. 3 ' 1415926 ' 
in2=S^ = ^ =.6931472, 

i=l 1 

<» r_i vo* 

In 3=S V = 1 .09861 2. 

These sequences are denoted by .1,, A,, and A z 
respectively. 

The first two examples are alternating rational 
sequences and tables (> and 7 show the results for Q 
and 8 X to be as predicted by theorems 3 and 9. 
1>(A), r D(A 2 ) J £ 2 (A l ) 9 and 6 2 {A 2 ) converge much 
more rapidly and apparently T) accelerates the 
convergence better than E 2 by a factor in \jn. 

Table 6 



II 


.1, 


CUi) 


fi(.l,) 


'D(.li) 


ft(il,) 


1 

2 
3 
4 
5 

6 


l.o 

2. 6666667 
3.4666667 
2.895238J 

3. 3396825 

2.9760462 
3. 2837385 

3. 01 707 is 
3. 2523H5!) 
3.0418396 

.099753 


3. 1466666 
3. L40136 
3. i421516 
3.1413354 
3.1417266 

3. 1415163 
3 1 U6392 


3. 1666667 
3. L333333 
3. 1452381 
3. L396825 
3.1427129 

3. L408813 

3. 11207 is 
3. 1412548 


3. 1417367 
3. L415567 

3. 1116012 

3. 1 415882 
3. 1415946 


3. 1423423 
3. L413919 
3. L416627 
3. 1 116634 

3. 11 looor. 

3. 1 U58! I 


7 






8 






9 






10 










error of last 
term 


.0000466 


. 0003378 


.01100020 


0000072 



Tabi e 7 



n 


A 2 


CUi) 


SiUi) 


^>( 1 


e*(Ai) 


1 
2 
3 
4 
5 


1.0 
.5 
. 83333333 

. 58333333 
. 78333333 


. 69444444 
.69270833 

. 60333333 
. 09305555 
. 69319726 


.7 

.60017010 
00144444 
. 69242424 
. 69358974 


.69318182 
.09313724 
.09315068 
.0! )3 14 573 
. 69314783 


. 69333333 
.69308943 
.69316940 
.69313724 
.69315212 


6 

7 
8 
9 
10 


. 61666666 
. 75952379 
. 63452379 
. 7456)3490 
. 64563490 


.69311755 
. 69316576 
. 69313490 
. 69315556 
. 69314123 


. 69285712 
. 69334732 

.69300331 
. 69325395 
. 69306573 


.69314682 
.09314734 
.69311700 
.00314723 
.69314713 


.69314448 
.69314872 
.69314621 
.6931 1778 
. 69314677 


11 
12 
13 
14 


. 73654400 
.65321067 
. 73013375 
.65870517 
72537184 


. 69315149 
. 69314395 
.69314961 
.69314529 
69314865 


.69321007 
.69300070 
.69318793 
.69311377 
69317488 


.69314719 
.69314715 

.09314719 


. 69314745 
.69314698 
.003117:51 

.6931 1707 


15 






16 


. 6(i287184 
.72169538 
. 66613982 




.69312394 






17 








18 










error of last 
term 


.0220074 


.0000015 


.0000232 


.00000001 


.00000010 



It is remarkable that -k can be found so accurately 
from a few terms of the simple sequence A x . Shanks 
considers £\{A^) and tt is given to 8 places by 
the first term of £\(Ai). It would require ap- 
proximately 40,000,000 terms of A x to give the 
same accuracy. 
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The third example is a "mixed" sequence. All 
of the transformed sequences diverge, but they are 
semiconvergent or asymptotic [3, pp. 520, 536] se- 
quences for In 3. In this example £ ana " ^ show 
a very marked improvement over £ x and £ 2 . It 
appears that round off has begun to affect T)(A 2 ). 
Although r D{A 2 ) was computed using double pre- 
cision, A 2 itself was computed to only eight signif- 
icant figures. 

Table 8 



n 


A, 


CUz) 


ffiMa) 


7>Uj) 


fiUi) 


1 
2 
3 
4 
5 

6 

7 

8 



*2.0 



2. 6666667 
-1.3333333 

5. 0666667 
-5.6 
12. 685714 
-19.314286 
37. 574603 

-64. 825397 
121. 35642 

120. 254 
. 801388 


1.0 

l.iiiini 

1.0952381 
* 1.0986667 
1.1008547 

1 . 0938776 
1. 1067668 
1 . 0852653 
1. 1202030 


1.0 

1.1428571 
1 . 0066666 
* 1.1282052 
1.0666666 

1.1368423 

1.0493503 
1.1657143 
1.0031748 
1.2391205 


1.0980392 
1.0980392 
1.1047618 
1.0991870 
1.0984127 

1.0986935 
* 1.0985975 


1.0909091 

1. 1014493 
1.0970464 
1.0997246 
1. 0976744 

1. 0995076 

* 1.0976739 

1.0996695 


9 






10 






11 








12 










error of last 

term 

*min error... 


.0215908 
. 0000545 


. 140508 
. 029593 


.0000147 
. 0000147 


. 0010573 
. 0008954 



The sequence A x may be modified by repeating 
each term in the sequence twice. Call the resulting 
sequence A[. The application of C> T>> &i an d 
£ 2 to A[ give some unusual results. It is seen 
that 



C{A[)-- 



--A[ 



[a n \a n = °°, 



1,2, 



8 2 (A[) is of the same type as A[ with each term 
repeated twice. The odd terms of T){A[) agree with 
those of £ 2 {A[), but the even terms do not. The 
even terms of £ (A[) converge much more rapidly 
than the odd ones, by a factor of approximately 
1/n. Both £ 2 (A[) and T)(A[) converge more slowly 
than £ 2 (Ai) and 7) (^0. 

5.3. Three Related Sequences 
The following three sequences are considered: 
P 3 ={(.98)*+(.95) w }, 
P*={(lM) n +(.95) n }, 
P 7 ={ (1.96)*+ (1.9)"}. 



The purpose is to investigate the following two 
questions. Is it advantageous to increase the largest 
term in an exponential sequence before applying 
C? Is it advantageous to keep the larger terms 
away from 1 at the expense of increasing them? 
The answer to the first question appears to be some- 
times and the answer to the second appears to be 
no. 

The limit or antilimit of each of these sequences 
is zero. The first converges quite slowly and the 
other two diverge (see table 9). 

Table 9. Behavior of P-s, P$, and P7 



n 


P3 


Pe 


Pi 


10 

25 
50 

75 


1.416 

.8808 
.4411 
.2411 


837.3 
2.025X107 
4.1X1015 


1450 
2.955X107 
4.967X10'5 









Ps, C(P 3 ), C(Pz\ and C 3 (P 3 ) all converge 
slowly. There is an improvement by a factor of 
approximately 2 in the range considered for each 
further application of C- 

Theorem 2 states that CiPz) and C(P&) converge 
at the same rate while C(P 7 ) diverges. This is 
borne out by table 10. Q 2 and C 3 produce entirely 
different effects for each sequence. £ 2 (P 6 ) con- 
verges very rapidly to zero and C 3 (P 6 ) converges 
even more rapidly. Errors due to the loss of sig- 
nificant figures have affected the 24th term of 
C 2 (P 8 ) and the 14th term of C 3 (P 6 ). All calcu- 
lations here were made with 20 decimal digits. 

In P 7 , C(Pi), C 2 (Pi), and C 8 (P 7 ) the first term is 
the closest to zero. However it is apparent that the 
sequence formed by the nth terms of C k (Pi), 
k=0, 1, . . . converges quite rapidly to the antilimit 
of P 7 . 

5.4. Effect of Lower-Order Terms 

To illustrate the effect of lower order terms we 
consider 

P 8 ={2*+(1.2) w }, 

which may be compared to P x in section 6.2. 



Table 10. Effect of C k on the sequences P 3 ={(.98)»+(.95)»|, Pg={ (1.96)«-h(.95)»|, P 7 ={(1.96)»+(1.9)« 



n 


C(Pz) 


C 2 (Ps) 


C 3 (Ps) 


C(Pe) 


cm 


C 3 (Pe) 


Ci.Pi) 


CKPi) 


CKPi) 


1 
5 

10 
15 
20 

25 
30 
35 
40 
45 


.2991 
.2708 
.2384 
.2091 

.1827 

.1589 
.1376 
.1186 
.1016 
.0866 


-.0716 
-. 0684 
-. 0638 
-. 0584 
-. 0525 

-.0462 
-.0399 
-. 0337 
-. 0279 
-. 0227 


-.0873 
-.0923 
-. 1032 
-. 1269 
-.2017 

-4. 903 
. 1319 
.0490 
.0241 
.0131 


1.024 

.8350 
. 6462 
.5000 
.3869 

.2994 
.2316 
.1792 
.1387 
.1073 


-. 0671 
-. 0028 
-.57X10- 4 
-. 12X10-5 
-.25X10-7 

. 15X10-8 
. 15X10-7 
.65X10-7 
.94X10-5 
.75X10- 1 


-.55X10-3 
-. 13X10-5 
-. 75X10- 9 
-. 12X10-'° 
-.21X10-9 

-.19X10-8 
-.38X10-9 
-.98X10-6 
-. 12X10- 4 


.0059 

.0807 
2. 129 
55.83 
1456 

.38X105 
.98X106 
.25X10* 
.64X109 
.16X1011 


-.88x10-5 
-.12x10-3 

-.0031 
-.0769 
-1.901 

-46. 27 
-1109 

-.26X105 
-.61X106 
-.14X108 


.40X10-7 
.54X10-6 
. 13X10-4 
.32X10-3 
.0076 

.1726 
3. 852 
84.02 
1795 
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Tablf 11 


Effect c 


/ lower or 


lev terms 








// 


Pi 


Ps 


C(Pi) 


C(P%) 


r 2 (p.) 


cm\) 


CKPi) 


O(Ps) 


1 
to 

20 

?7 


4. 8426 
1018. 569885 
10,363,260.20 
1.326,449,998 


3. L626 
I0IS. 110177 
10,363,258.88 

1,326, 149,997 


2 S810 

4. 8408 
26. 826 
95. 645 


.8205 
4. 3071 
L'f,. 673 
95. 576 


2. 695 

1.282 
. 3477 


-.0719 

-.00i78 

-.00014 


2.711 
-.0072 
.0267 


2XHH 
3X10- 
5X10-" 













In (able 11 il is seen that the lower order terms 
disappear rather rapidly in P u but the 27th terms of 
C(P\) and C(Ps) on ly agree to two digits. Further 
C 2 {P\) and C 3 (Pi) are completely unrelated to £ 2 (P 8 ) 
and C 3 (Ps)- Thus it is seen that these lower order 
terms are significant even for CV*) and they affect 
the rate of convergence of C k (P), k^>l, drastically 
and adversely. 

5.5. Loss of Significant Figures 

It has been seen that terms witli a small effect on 
P can have a large effect on C k {P)> A related 
phenomenon is the loss of significant figures in the 
computation of C(P) and T>(P). If P is the con- 
stant sequence {p n \p n =l}, then C(P) and T)(P) 
may be written symbolically as 

i-q-D 
(i-i) 

Hence slight inaccuracies in the computation of 

CiP) and THJ.*) may distort the results badly if P 
is slowly varying. Often it is necessary to use double 
precision for tin 4 computation of '2)(/ > ) even though 
P is only given to six or eight digits. 

The following table shows where the loss of signifi- 
cant figures can occur. The first column of table 12 
defines the sequence P, the second column gives the 
approximate size of the terms of P in the range that 
the transformation is applied and the third column 
gives the transformation used. The next four 
columns give n, the size of the terms appearing in 



the numerator, the actual value of the numerator 

and the number of significant figures lost from P. 

Table 12. Loss of significant figures 



Sequence 


Size of 
terms 


Trans- 
forma- 
tion 


n 


Size of 

numei leal 

terms 


Numer- 
ator 


Signif- 
icant 

figures 
lost 


(1.96)»+(.95)» 


.1X108 

.4X10-' 

.5X10-3 

2 

1 




25 

30 
75 
all 
15 


.01 
. 0002 
.2X10-9 
10 

1 


.43X10-* 

.3X10- 10 

.2X10- 21 

10-M 

io- il 


20 


(.9) "+(.8) "cos (An)... 
(.9) "+(.8) »cos (An)... 
(.999)»-f(.99)"- 


9 

is 

19 


(.98)»+(.6)M-(.4)»-- 


11 
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